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Abstract 

A Higgs particle produced in association with a Z boson and decaying into weak 
boson pairs is searched for in 336.4 pb^^ of data collected by the L3 experiment 
at LEP at centre-of-mass energies from 200 to 209 GeV. Limits on the branching 
fraction of the Higgs boson decay into two weak bosons as a function of the Higgs 
mass are derived. These results are combined with the L3 search for a Higgs boson 
decaying to photon pairs. A Higgs produced with a Standard Model e"'"e~ Zh 
cross section and decaying only into electroweak boson pairs is excluded at 95% CL 
for a mass below 107 GeV. 



Submitted to Phys. Lett. B 



1 Introduction 



In the Standard Model, the Higgs mechanism generates the masses of elementary particles and 
stabihses the high energy behaviour of the electroweak interactions. However, the Standard 
Model does not predict the mass of the Higgs boson. Searches to date have focused on b 
quark decays of the Higgs, but the model predicts an increased branching fraction to massive 
vector boson pairs for a heavier Higgs. In some extensions of the Standard Model which predict 
multiple Higgs bosons, the lightest Higgs boson couples primarily to bosons, not fermions [1]. 
Results excluding these "fermiophobic" models were reported for Higgs decays into two photons, 
but for Higgs masses above 90 GeV, the decay into massive vector boson pairs dominates [2]. 

Several possible models predict the presence of a fermiophobic Higgs. As a point of reference, 
a simple extension of the Standard Model is chosen where the production cross section for 
e+e^ Zh, the Higgsstrahlung process, is kept at the Standard Model value. All direct decays 
into fermions are removed, and the resulting branching ratios favour decays into two photons for 
Higgs bosons with masses below 90 GeV, while massive vector boson pairs become important 
above 90 GcV. This set of assumptions is called the "benchmark fermiophobic scenario". The 
branching ratios to 77, W+W^, and ZZ predicted in this scenario are plotted in Figure 1. 

This Letter presents the search for a Higgs boson produced in association with a Z boson 
through the process e+e" — > Zh, followed by the decay of the Higgs to either a pair of W or 
Z bosons. Throughout the analysis we assume that the ratio between the branching fraction 
h — > WW* and the branching fraction h — > ZZ* is given by the Standard Model expectations. 
This assumption is valid for scenarios beyond the Standard Model satisfying the constraint 
p ~ 1 [3]. Consequently, the branching fraction h — > WW* is expected to be the dominant 
weak boson decay for all Higgs masses. Accordingly, the analyses focus on h — > WW* but 
the efficiency of the analyses for a h — > ZZ* signal is considered where appropriate. Since the 
Higgsstrahlung mass reach of the LEP data is less than 160 GeV, at least one of the weak 
bosons produced from the Higgs decay must be off its mass shell. The analyses assume the 
presence of one vector boson with a mass within ±10 GeV of its nominal mass and one with 
a much smaller mass. This topology is expected to be present in approximately 70% of signal 
events. 

2 Data and Monte Carlo Samples 

We analyse data collected with the L3 detector [4] at ccntrc-of-mass energies ^/s = 200 — 
209 GeV, for a total integrated luminosity of 336.4 pb^^. The data are grouped into six samples 
whose average centre-of-mass energies and corresponding integrated luminosities are listed in 
Table 1. 

The Higgs signal cross section and fermiophobic branching ratios are calculated using the 
HZHA Monte Carlo generator [5]. For efficiency studies, Higgs events are generated using 
PYTHIA [6] for Higgs masses between 80 GeV and 120 GeV, with a mass step of 5 GeV up 
to 100 GeV and 1 GeV thereafter. The e+e" ZZ, e+e" Ze+e", and e+e" qq{j) 
processes are simulated with the PYTHIA generator. The KK2f [7] generator is also used for 
the e+e" — > qq(7) process. The KORALW [8] generator is used for the e"'"e~ — > W''"W~ process 
except for the e''"e~ — > e~Peqq' final state, which is simulated using EXCALIBUR [9]. Hadron 
production in two-photon interactions is modelled using the PHOJET program [10]. 

The L3 detector response is simulated using the GEANT program [11], which takes into 
account the effects of energy loss, multiple scattering and showering in the detector. The 
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GHEISHA package [12] is used to simulate hadronic interactions in the detector. Time depen- 
dent detector inefficiencies, as monitored during tfie data acquisition period, are also simulated. 

3 Analyses 

Of the possible final states from the e+e^ — > Zh — > ZWW* — > 6 fermions process, this Letter 
describes the results of analyses for six of them, covering 93% of the total branching ratio. 
The six decay channels searched are: Zh qqqqqq, Zh — > qqqq^i/, Zh — > qqii/iu, Zh — > 
i/i/qqqq, Zh uuqqiu, and Zh — > ££qqqq ^\ The analyses adopt a common procedure of 
preselection, followed by a selection using a neural network approach and the production of a 
final discriminant variable. 

The first step is to apply a set of preselection cuts to remove the background events most 
different from the signal. Important preselection cuts are made on visible energy and counts 
of tracks and clusters. The number of neutrinos in a channel determines the visible energy 
window, while the number of jets provides a basis for cuts on the number of charged tracks and 
calorimeter clusters. For channels containing leptons, important lepton identification cuts are 
applied. Two-photon events and e"'"e~ —>■ qq(7) events containing one or more photons from 
initial state radiation are suppressed using cuts on the fraction of deposited in detectors 
near the beam line and by cutting events where the missing momentum vector points near the 
beam pipe. These cuts also reduce the number of e+e" ei/qq background events accepted 
by the analysis. 

The background processes which are more similar to the signal are removed using one or 
more neural networks [13]. Where possible, preselected events are subjected to a constrained fit 
to the expected event topology before constructing the neural networks. Variables considered 
in the networks include the energies of the most and least energetic jets, the minimum angle 
between any two jets, the minimum angle between any jet and any lepton, and the reconstructed 
masses of any Z or W bosons identified in the event. 

After applying cuts on the neural networks, the analyses produce final distributions of the 
selected signal, background, and data, generally using a discriminant combination of the neural 
network outputs and a reconstructed Higgs mass. 

In the following sections wc outline the selection procedures for the six channels. Complete 
details are available in Reference 14. 

3.1 Zh — > qqqqqq 

The qqqqqq analysis searches for the case when the W and Z bosons decay into hadrons. The 
primary backgrounds to this channel are ZZ qqqq, W"'"W~ —>■ qqqq and e"'"e~ — qq(7) 
events. This six-jet signature is also produced by the Zh — > ZZZ* — > qqqqqq process, which 
represents 34% of the total h ZZ* signal. The efficiency for the h ZZ* signal is included 
in the analysis, effectively adding 15% to the expected rate relative to using only h WW*. 

Events with full energy and large hadronic content are selected. The analysis uses a fit of 
the event to a six-jet topology using the Durham algorithm [15]. Events with a poor match to 
the six-jet hypothesis are rejected. 

^^This simplified notation, covering all possible combinations of flavour and charge, is adopted throughout 
this Letter. 
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After preselection, each six-jet event is subjected to a kinematic fit, which requires momen- 
tum and energy conservation. Following the fit, the pair of jets with invariant mass closest to 
rriz is chosen as the Z candidate. Of the remaining four jets, the dijet pair with invariant mass 
closest to mw is assigned to be the W candidate and the remaining pair is identified as the W*. 

The analysis uses three neural networks, one for each of the e+e~ — > W+W~, e+e" — > ZZ, 
and the e"^e~ — > qq(7) backgrounds. All three neural networks use the same eleven inputs and 
have twenty hidden nodes and one output node. The input variables are chosen according to 
the specific features which differentiate signal from background: 

• Typical QCD features in the background are: unequal repartition of the energy among 
the six jets, jets with low multiplicity, small invariant masses and small angles between 
closer jets. The following six input variables are therefore chosen: the energies of the 
most and least energetic jets, Ejl^^ and EJH^, the minimum number of tracks in a jet, 
njgj", the minimum angle between jets, OJ^"'' and the values of the Durham algorithm 
parameter logy45 and logyse- The value ymn is the threshold for which the DURHAM 
algorithm reconstructs a m-jet event as a n-jet one. 

• The event must be kinematically inconsistent with the production of only two on-shell 
boson pairs. This leads to the choice of the following five input variables: the chi-square 
of a kinematic fit to equal-mass boson pair-production, Xww^ mass obtained in the fit, 
mlg, the masses of the Z and W candidates after kinematic fit to the signal hypothesis, 
m^** and m{^, and the angle between the decay planes of the W and W* candidates, 
o^ww* • 

As an example of the input variables to the neural networks. Figure 2a presents the dis- 
tribution of for data, background and signal Monte Carlo. Table 2 gives the numbers of 
signal and background events expected and data observed after cuts on the neural network out- 
put. The final variable produced by the analysis is a discriminant combining the three network 
outputs and the reconstructed Higgs mass after the kinematic fit and Z candidate assignment. 
Distributions of this final variable for data, background and signal are plotted in Figure 3a. 

3.2 Zh ^ qqqq£i/ 

In the qqqq^i/ channel, the Z decays into hadrons, while one W decays into hadrons and the other 
decays into leptons. The different lepton flavours naturally define three different subchannels: 
qqqqez/, qqqq^z/, and qqqqrz/. Further, the difference between leptons coming from the W and 
from the W* doubles the number of subchannels. In one set of signatures, the W decays into 
hadrons and the W* decays into leptons, which means the lepton and neutrino energies are 
small, as is the missing energy of the event. In the other set, the W decays into leptons and the 
W* decays into hadrons, yielding a high energy lepton and large missing energy. The major 
backgrounds are W pair production (especially e"'"e~ W"'"W~ — > qqri^), Z pair production, 
and the c+c^ fiq(7) process. 

Events are classified into a subchannel according to the most energetic lepton in the event. 
For the qqqqei^ and qqqq//i/ cases, the subchannels are separated using the ratio of the lepton 
energy to the visible energy. In the subchannels where the lepton-neutrino pair comes from the 
on-shell W, both the lepton and neutrino energies should be large. For the subchannels where 
the lepton and neutrino come from the W*, the converse is true. In the qqqqrz/ channels, the 
initial lepton energy is difficult to reconstruct, so the subchannels are separated using only the 
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visible energy. The candidate lepton in the event is required to meet minimum identification 
criteria to ensure isolated and well reconstructed leptons. 

The analysis selects events with significant hadronic content and an appropriate amount of 
visible energy depending on the subchannel. A major source of background for this channel is 
the W+W~ qqlu process. To reduce it, an event is rejected if its hadronic mass is less than 
90 GeV. After preselection, the major remaining background is e+e~ — > W+W~, particularly 
e''"e~ W+W~ qqqq events with a lepton from hadron decays. 

Two neural networks with ten input nodes are prepared for each lepton fiavour to remove 
the e+e" W+W^ and e"'"e~ euqq backgrounds. Besides the common input variables E™t^ , 
^jet \ ^Ti^^ Xww' ''^z** ^^"^ log 1/34) the networks differ in the hypothesis used in associating 
the lepton-neutrino or quark- antiquark pairs to the W or the W*^*) bosons. For each network, 
four additional variables, constructed after a kinematic fit to the qqqq^i^ topology, are used: 
the invariant masses of the lepton-neutrino and quark- antiquark candidate systems, m£* and 
m^q , as well as their respective momenta, p^** and p^** . These momenta are expected to be 
small for the signal hypothesis. Twenty hidden nodes and one output node are used for each 
network. The distributions in data, background and signal Monte Carlo of the variable m^r^* 
used in the neural networks are presented in Figures 2b and 2c for the hypotheses of leptonic 
decays of the W and W*^*-' bosons, respectively. The mimbers of events expected and observed 
in this channel after cuts on the neural network outputs arc listed in Table 3. 

The mass of the Higgs boson is reconstructed from the recoil against the two jets which 
form the Z candidate. The final result of the analysis is a discriminant variable combining 
the output of the neural network and the reconstructed mass separately for each subchannel. 
Figure 3b shows the distributions of the final variable with all six subchannels combined. 

3.3 Zh ^ qq^i/^i/ 

In the qqlvlv channel, the Z decays into hadrons while both W bosons decay into lepton- 
neutrino pairs. The analysis requires two identified leptons, one with more than 12 GeV of 
energy and the second with energy greater than 10 GeV. The two neutrinos in the signal 
signature imply a visible energy window between 50% and 85% of \/s. The analysis starts from 
events with hadronic content and rejects the two-photon and e''"e~ — > qq(7) backgrounds as 
described above. 

The dominant backgrounds after preselection are the c+c^ W+W~ and e+e~ — > ei/qq 
processes. They are suppressed by using a single neural network based on seven kinematic in- 
put variables, with sixteen hidden nodes and one output node. No constrained fit is performed 
for the Zh — > qqlvlv hypothesis because the two neutrinos render the technique ineffective, 
but some of the input variables are similar to the ones employed in the previous sections: 
log(y23/?/34)) the hadronic invariant mass, rrihad, which shoTiId be consistent with rriz, the min- 
imum angles between each of the leptons and the closest jet, expected to be small in the 
case of quark decays, the invariant masses of the lepton-neutrino system and of the remaining 
of the event after a kinematic fit to the W+W" — > qqiv hypothesis, m{** and m{.ll^, and the 
angle between leptons, 9ii, which should be high in the background case. Distributions of nihad 
in data, background and signal Monte Carlo are presented in Figure 2d. The numbers of events 
expected and observed in this channel after a cut on the neural network output are listed in 
Table 4. 

After selection, the events are divided into two groups: events where neither lepton is 
identified as a tau and events where at least one of the leptons is identified as a tau. Most 
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background events fall into the second group. Final discriminant variables are constructed for 

these two subchannels separately, using the reconstructed Higgs mass and the neural network 
output. The Higgs mass is reconstructed by scaling the jet masses and energies by a common 
factor until the dijet mass is equal to mz, and then calculating the recoil mass against the dijet. 
Figure 3c shows the distributions of the final variable for data, background and signal, with 
the two groups of events combined. 

3.4 Zh '''i^qqqq 

In the case where the Z decays into two neutrinos and both W's decay into hadrons, the 
signature is four jets plus missing energy. The missing mass should be the Z mass and the 
visible mass the Higgs mass. As in the qqqqqq channel, the Zh ZZZ* — > z/z/qqqq events 
are included along with the h — > WW* signal. The selection accepts events where either the 
radiated Higgsstrahlung Z or the Z from the Higgs decays into neutrinos, but not the Z* — > uu 
case because of insufficient missing energy. The accepted signatures comprise 20% of the total 
Zh — > ZZZ* branching fraction, and their inclusion increases the expected signal rate by 15%. 
Most background comes from Z and W pair production (especially e"'"e~ — > W'*'W~ — > qqri/), 
and the e^e~ — qq(7) process. 

The analysis starts from events with significant hadronic content as well as substantial 
missing energy, suppressing e'*"e~ — > qq(7) and two-photon processes. Events are constrained 
into a four-jet topology using the Durham algorithm. Events where any of the four jets contains 
no charged track or has an energy less than 6 GeV are rejected. 

The preselected events are subjected to a constrained kinematic fit assuming a balanced 
event with four jets and an invisible Z. The two jets with the invariant mass closest to mw 
are considered as the W candidate dijet. Three separate neural networks are created and used 
to remove background processes. These three networks, based on ten input variables, twenty 
hidden nodes and one output node, are trained against Z and W pair production and the 
e"'"e~ qq(7) process. The eight input variables used in the networks are: the recoil mass of 
the event, rrirec, which should be consistent with mz for signal events, E™^^, EJ^]^, 0^^"", rriyf, 
log 1/23, log 1/34 and a^w*- The angle a^w* corresponds to the angle aww* calculated in the 
rest frame of the hadronic system. Distributions in data, background and signal Monte Carlo 
of a^y^t are presented in Figure 2e. Table 5 gives the numbers of signal and background events 
expected and data observed after cuts on the neural network outputs. 

The final variable is constructed by combining the outputs of the three neural networks 
with the reconstructed Higgs mass in a single discriminant. Figure 3d shows the distributions 
of this final variable for data, background and signal. 

3.5 Zh — > i/i/qq£i/ 

In the i/i/qqiiy channel, the Z decays into neutrinos and the WW* pair decays into two quarks, 
a lepton and a neutrino. The major backgrounds for this process depend on lepton flavour and 
include W and Z pair production, ez/qq and Zee final states, as well as hadronic two-photon 
processes. As in the qqqqfz/ channel, the signal divides into six subchannels as a function of 
the lepton flavour and origin. The same variables and lepton identification requirements are 
used to separate the subchannels. 

The small visible energy in this channel makes the suppression of two-photon processes 
particularly important. Stringent requirements against this background are applied and events 
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with moderate hadronic content are selected. 

After preselection, the most important background is e+e" — > W+W" qqiu (particularly 
qqrz/). This background is rejected using several variables including the average reconstructed 
W pair mass, m^^. This variable is constructed by scaling the energies and masses of the two 
jets by a common factor so that the sum of the jet energy is y/s/2, as would be the case in 
a e+e~ — > W+W" — > qqiu event. Using these rescaled jets, the W mass is calculated as the 
average of the dijet invariant mass and the invariant mass of the lepton and the missing energy 
vector. 

Each subchannel uses one neural network with eight input variables, twenty hidden nodes 
and one output node in order to remove the W+W~ and ei/qq backgrounds. The neural 
network input variables include m^^, the dijet invariant mass and the missing invariant mass. 
In addition, most of the remaining kinematic information in the event is used: the lepton energy, 
the hadronic energy, the minimum angle between the lepton and jets, the angle between the 
lepton and the dijet plane, and the angle between lepton and missing energy vector. The 
numbers of events predicted and observed after cuts on the neural network outputs are listed 
in Table 6. 

There are insufficient constraints to fully reconstruct the Higgs mass, so the visible mass is 
used as the final variable after a cut on the neural network. Figure 3e shows the distributions 
of the final variable with all six subchannels combined. 

3.6 Zh ^ iiqqqq 

In this channel, the Z decays into a pair of electrons or muons and hadronic W decays are 
considered. The event signature is two energetic leptons with the invariant mass of the Z, two 
energetic jets with an invariant mass near mw and two lower energy jets. This signature is 
also produced by 6.5% of the h — > ZZ* signal events. 

The events are separated into two subchannels according to the lepton flavour. The analysis 
selects balanced energy events with significant hadronic content. 

The preselected events are fit to a topology with two leptons and four jets and subjected 
to a kinematic fit requiring energy and momentum conservation. The eeqqqq subchannel uses 
a neural network which is trained against the e'^e" — > ZZ and e"'"e~ — > Zee backgrounds. 
The ///xqqqq subchannel uses two neural networks, one trained to discriminate against the 
e"'"e~ ZZ process and the other against e"'"e~ — > W^W~. The neural networks have six 
input nodes, twenty five hidden nodes and one output variable. The six input variables are: 
Ejl^^, E™j^, OJ-j^, ^mm^ logy34 and the invariant mass of the dilepton system, which should be 
consistent with mz- The number of background and signal events expected and data events 
observed after cuts on the neural network outputs are given in Table 7. 

The final variable is a discriminant which combines the output of the ZZ rejection neural 
network with the Higgs mass reconstructed from the recoil against the lepton pair. Figure 3f 
shows the final variable distributions with the two subchannels combined. 

4 Results 

The sensitivity to a Higgs boson decaying into weak boson pairs is enhanced when all the pre- 
vious analyses are combined. This combination considers each final variable as an independent 
Poisson counting experiment. A single statistical estimator from all channels is then built. This 
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statistical estimator is the log- likelihood ratio [16]: 



21ng = 2^Sj - 2^njln ( 1 + ^ 




where rii is the number of data events observed in the i-th bin of the final variable, Sj is the 
number of expected signal events and bi is the number of expected background events. The 
values of n^, 6j and Sj depend on the Higgs mass hypothesis. The —2liaQ values expected in 
the presence of background alone or both signal and background are determined by replacing 
rii by bi or Si + bi, respectively. The plot of the —2\nQ as a function of the mass hypothesis 
is shown in Figure 4. The dark band surrounding the data gives the size of the systematic 
uncertainties. The systematic dependence of the —2lnQ distribution on detector quantities 
such as calorimeter energy scale and tracking efficiency is determined by shifting correlated 
variables in the analysis and propagating the effects to the signal and background hypothesis 
distributions. Other factors, such as signal and background Monte Carlo statistics and cross 
section uncertainties, are also included. Table 8 details the different sources of systematic and 
their contribution to the total systematic uncertainty on —2 In Q. 

No significant excess indicating the presence of a Higgs boson decaying into WW* or ZZ* 
is observed in the data, which would manifest as a significant dip in the —2 In Q distribution. 
Confidence levels for the absence of a signal are hence derived [16]. Figure 5 shows the observed 
and expected 95% confidence level (CL) limits as a function of the Higgs mass, assuming the 
Standard Model production cross section. In the assumption of BR(h WW*)+BR(h — >■ 
ZZ*) = 100%, a Higgs boson with mass less than 108.1 GeV is excluded at 95% CL. Assuming 
the value of BR(h — > WW*)+BR(h — > ZZ*) of the benchmark fermiophobic scenario, calculated 
with the HDECAY program [17], the observed exclusion region is 83.7 GeV<mh< 104.6 GeV 
with a region between 88.9 GeV<mh< 89.4 GeV which can be excluded only at 93% CL. 

Model independent fermiophobic results can be derived by scanning the relative branch- 
ing fractions of h — > 77 and h WW*. The branching fractions into boson pairs can be 
conveniently parameterised in the form: 



R^j represents the fraction of fermiophobic decays into photon pairs, and ranges from zero to 
one, while BRbosons represents the total Higgs branching fraction to pairs of gauge bosons. The 
scan combines these h WW* results with the previously published L3 h ^ 77 results [18], 
determining the 95% CL exclusion for BRbosons at each point in the mh versus R.yy plane. The 
full scan results are presented in Figure 6. In the benchmark scenario, the fermiophobic Higgs is 
limited at 95% CL to have > 108.3 GeV, compared to an expected hmit of mh > 110.7 GeV. 
These results represent a significant extension of the > 105.4 GeV limit obtained in the 
photonic channel [18]. These results also exclude at 95% CL any model with < 107 GeV for 
BRbosons — 100% and any value of R-y-y, assuming the Standard Model production cross section. 



BRbosons = BR(h^77) + BR(h^WW*) + BR(h^ZZ*), 
R^^ = BR(h 77) / BRbosons- 
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(GcV) 


199.5 201.8 203.1 205.0 206.5 208.0 


Luminosity (pb~^) 


82.8 37.0 8.8 68.9 130.4 8.5 



Table 1: Average centre-of-mass energies and corresponding integrated luminosities. 



Zh (1(1(1(1(1(1 


Data 


Background 


\\'+\\'- 


ZZ 


cici(:.0 


Signal 


Preselection 


1886 


1870.1 


1274.4 


104.6 


488.9 


16.6 


Selection 


443 


446.0 


347.7 


44.1 


54.0 


14.4 



Table 2: Number of events observed in data by the Zh qqqqqq analysis, compared with 
the Standard Model expectations. The Monte Carlo breakdown in different processes is given. 
Signal expectations are given for mh = 105 GeV in the fermiophobic benchmark scenario. 



Zh — qqqqez/ 


Data 


Background 


w+w- 


ZZ 


qq(7) 


Signal 


Preselection 


13 


12.5 


4.1 


2.3 


2.8 


2.3 


Selection 


6 


5.3 


1.5 


1.4 


1.1 


2.0 


Zh —>■ qqqq/iz/ 


Data 


Background 


w+w- 


ZZ 


qq(7) 


Signal 


Preselection 


3 


7.3 


5.0 


1.4 


0.9 


1.7 


Selection 


1 


3.0 


1.7 


0.8 


0.5 


1.5 


Zh qqqqrz/ 


Data 


Background 


w+w- 


ZZ 


qq(7) 


Signal 


Preselection 


251 


202.4 


137.5 


13.9 


43.4 


1.6 


Selection 


41 


41.2 


24.4 


2.8 


8.6 


1.3 



Table 3: Number of events observed in data by the Zh qqqqez/ analysis, compared with 
the Standard Model expectations. The Monte Carlo breakdown in different processes is given. 
Signal expectations are given for = 105 GeV in the fermiophobic benchmark scenario. 
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Zh — > qqiuii' 


Data 


Background 


W+W" 


ZZ 


ez/qq 


Signal 


Preselection 


23 


21.4 


10.7 


1.6 


7.6 


0.7 


Selection 


2 


1.4 


0.5 


0.5 


0.2 


0.5 


Zh —f i[{[^Tu(u 


Data 


Background 


\\+\\- 


ZZ 




Signal 


Preselection 


96 


91.7 


57.3 


6.4 


20.8 


0.6 


Selection 


28 


23.1 


13.1 


3.5 


4.6 


0.4 



Table 4: Number of events observed in data by the Zh — > qqiviv analysis, compared with 
the Standard Model expectations. The Monte Carlo breakdown in different processes is given. 
Signal expectations are given for mh = 105 GeV in the fermiophobic benchmark scenario. 



Zh — > z/z/qqqq 


Data 


Background 


w+w- 


ZZ 


ez/qq 


qq(7) 


Signal 


Preselection 


451 


439.7 


263.2 


25.3 


73.6 


75.9 


4.0 


Selection 


41 


47.9 


41.5 


4.9 


4.4 


6.9 


3.3 



Table 5: Number of events observed in data by the Zh z/z/qqqq analysis, compared with 
the Standard Model expectations. The Monte Carlo breakdown in different processes is given. 
Signal expectations are given for rrih = 105 GeV in the fermiophobic benchmark scenario. 



Zh — > pi/qqiv 


Data 


Background 


w+w- 


ZZ 


euqq 


Signal 


Preselection 


21 


11.3 


5.7 


1.2 


3.9 


1.0 


Selection 


2 


1.7 


1.0 




0.7 


0.8 


Zh uiyqq{£iy)* 


Data 


Background 


w+w- 


ZZ 


ei^qq 


Signal 


Preselection 


39 


44.4 


34.4 


3.6 


5.4 


0.6 


Selection 


7 


5.3 


3.9 


0.5 


0.9 


0.6 



Table 6: Number of events observed in data by the Zh — > i/i/qqiu analysis, compared with 
the Standard Model expectations. The Monte Carlo breakdown in different processes is given. 
Signal expectations are given for mh = 105 GeV in the fermiophobic benchmark scenario. 



Zh — eeqqqq 


Data 


Background 


Zee 


ZZ 


Signal 


Preselection 


4 


6.8 


1.3 


5.3 


0.4 


Selection 


2 


3.0 


0.7 


2.3 


0.4 


Zh — > /i/xqqqq 


Data 


Background 


w+w- 


ZZ 


Signal 


Preselection 


13 


13.7 


7.2 


5.5 


0.4 


Selection 


3 


1.2 


1.1 


2.7 


0.1 



Table 7: Number of events observed in data by the Zh ££qqqq analysis, compared with 
the Standard Model expectations. The Monte Carlo breakdown in different processes is given. 
Signal expectations are given for mh = 105 GeV in the fermiophobic benchmark scenario. 



14 





Analysis Zh 


Source of Systematics 


qqqqqq 


qqqqi?z/ 




z/i/qqqq 




iiqqqq 


Tracking effciency 


2.1% 


1.7% 


0.4% 


0.2% 


1.8% 


2.1% 


Energy scale 


1.3% 


2.9% 


2.0% 


6.0% 


5.3% 


2.2% 


Lepton identification 




3.8% 


2.8% 




4.4% 


1.8% 


Event shape 


1.2% 






2.6% 






Signal Monte Carlo statistics 


2.0% 


3.0% 


4.0% 


2.0% 


2.0% 


2.0% 


Background Monte Carlo statistics 


0.1% 


1.2% 


0.1% 


0.1% 


1.1% 


0.1% 


Background cross sections 


0.4% 


1.2% 


0.8% 


0.7% 


1.1% 


0.4% 


Total 


3.4% 


6.1% 


5.3% 


6.9% 


7.6% 


4.1% 



Table 8: Systematic uncertainties on the —2\nQ distributions for the different channels. 




m, (GeV) 

Figure 1: Higgs branching fractions predicted in the benchmark fermiophobic scenario. 
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Figure 2: Distributions of variables used in input to the neural networks for data, background 
and signal Monte Carlo at = 105 GeV. All centre-of-mass energies and subchannels are 
combined, b) refers to leptons originating from a W boson while c) refers to leptons originating 
from a W* boson. Signal Monte Carlo is magnified by the factor indicated in the figures. 
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Figure 3: Distributions of the final variable for data, background Monte Carlo and signal Monte 
Carlo at mh = 105 GeV. All centre-of-mass energies and subchannels are combined. The final 
variable is displayed in bins of logio( background )- 
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Figure 4: Log-likelihood ratio plot for the combined h — > WW*, ZZ* search. The dashed lines 
represent the value of the expected background-only and signal+background distributions. The 
shaded regions around the background-only line indicate the ±lcr and ±2cr regions. The solid 
line indicates the observed values. The bands around the observed sensitivity represent the 
effects of systematic uncertainties. 
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m, (GeV) 

Figure 5: Excluded values at 95% CL of BR(h ^ WW*) + BR(h ^ ZZ*) as a function of the 
Higgs mass (solid line), in the assumption of the Standard Model production cross section. The 
expected 95% CL limit (dashed line) and the fermiophobic benchmark prediction (dotted line) 
are also presented. The Standard Model prediction for BR(h — > WW*) is 8% at nih = 115 GeV, 
falling below 1% for mh < 100 GeV. 
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m, (GeV) 

Figure 6: The 95% CL limit for BRbosons as a function of rrih and FLy^. The sohd hues indicate 
the borders between regions of exclusion. The crossing point between the BRbosons ~ 100% 
line and the dashed line provides the lower limit on the Higgs mass in the benchmark scenario: 
rrih > 108.3 GeV. The minimum value of mh on the BRbosons = 100% contour gives a lower 
mass limit for any model where the Higgs decays exclusively into electroweak boson pairs: 
mh > 107 GeV. 
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